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Abstract 

We discuss the fragmentation of a heavy quark to a baryon containing two heavy 
quarks of mass mg ^> Aqcd- In this limit the heavy quarks first combine perturbatively 
into a compact diquark with a radius small compared to 1/Aqcd ? which interacts with 
the light hadronic degrees of freedom exactly as does a heavy antiquark. The subsequent 
evolution of this QQ diquark to a QQq baryon is identical to the fragmentation of a heavy 
antiquark to a meson. We apply this analysis to the production of baryons of the form 
ccq, bbq, and bcq. 
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The spectroscopy and interactions of baryons consisting of two heavy quarks and one 
light quark simplify in the limit that the heavy quark masses uiq tend to infinity. This 
is because the heavy quarks are bound into a diquark whose radius tqq is much smaller 
than the typical length scale 1/Aqcd of nonperturbative QCD interactions. In the limit 
r QQ 1/Aqcd the heavy diquark has interactions with the light quark and other light 
degrees of freedom which are identical to those of a heavy antiquark. Hence as far as these 
light degrees of freedom are concerned, the diquark is nothing more than the pointlike, 
static, colour antitriplet source of the confining colour field in which they are bound 

One immediate result of this limit is that the spectrum of such "doubly heavy" baryons 
is related to the spectrum of mesons containing a single heavy antiquark |l[ . It also follows 
that the form factors describing their semileptonic decays may be related to the Isgur-Wise 
function, which arises in the semileptonic decay of heavy mesons . In this note we will 
apply the same symmetries to the nonperturbative dynamics which governs the production 
of such states via fragmentation processes. We will use our results to estimate the produc- 
tion rates for baryons of the form ccq, bbq and bcq; however, we note that, especially in 
the cc system, the heavy diquarks are not particularly small relative to 1/Aqcd, so there 
may well be sizeable corrections to our results. 

The fragmentation of a heavy quark Q into a QQq (or QQ'q) baryon factorises into 
short-distance and long-distance contributions. The heavy quark first fragments into a 
heavy diquark via a process which is perturbatively calculable. In fact, the amplitude may 
be trivially related to that for the fragmentation of Q into quarkonium QQ. The subsequent 
fragmentation of the diquark QQ to a baryon is identical to the fragmentation of a Q to 
a meson Qq; this information may be obtained from experimental data on production of 
heavy mesons. 

We will begin with the case of baryons of the form QQq, in which the two heavy 
quarks have the same flavour. For concreteness we will discuss the production of baryons 
with two charm quarks ccq; the extension to bottom baryons bbq will be trivial. As the 
colour wavefunction of the charm quarks is antisymmetric and the quarks are taken to be 
in the ground state S'-wave, the spin wavefunction must be symmetric. Hence the cc can 
only form a spin-1 diquark, which we shall denote by (cc). The (cc) can then fragment 
either to a spin-| baryon, which we shall denote by E cc , or to a spin-| baryon, which we 
shall call S* c . 

The initial short distance fragmentation process c — > (cc)c is analogous to that for the 
fragmentation into charmonium, c — > ipc, which has been shown by Braaten, Cheung and 
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Yuan H to be calculable in QCD perturbation theory. The Feynman diagrams responsible 
for (cc) production are shown in fig. 1. The computation, which we shall not repeat, follows 
directly that outlined for c — > ijjc in [|]. In fact, after some rearrangement the form of the 
fragmentation function D c ^( cc )(z) is exactly the same as that for D c ^(z), except for an 
overall normalisation factor. Since the charm quarks in the (cc) are in an overall colour 
3 rather than a singlet, there is a colour factor of 2/3 instead of 4/3 in the amplitude. 
The colour wavefunction of the diquark is |(cc) a ) = \t a \, c |c 6 ) |c c ) where a, 6, c are colour 
indices. The factor of 1/2 in the normalisation of the diquark state is canceled by a factor 
of 2 in the matrix element from fermi statistics, giving an overall 2/3 in the amplitude. 
Squaring and summing over final colours then gives the desired fragmentation function, 
renormalised at the scale \i = 3m c : 

D c ^ {cc) (z,^ = 3m c ) = jL ^™^ ^ a s (3m c ) 2 F(z) , (1) 

where 

F( z ) = Z ^~^e (16 - 32z + 72z 2 - 32z 3 + bz 4 ) . (2) 

Here -R( cc )(0) is the nonrelativistic radial wavefunction at the origin for the perturbatively 
bound diquark. Unlike the case of charmonium, there is no physical "decay constant" 
/( cc ) to which it may be related. However we may naively scale i?( cc )(0) from charmonium 
by noting that in a hydrogen-like potential -R(O) is proportional to (CVa s ) 3 / 2 where the 
colour factor Cp is | for cc and | for cc. Hence we expect that |i?( cc )(0)| 2 ~ |-R^(0)| 2 /8 ~ 
(0.41 GeV) 3 . 

The complete fragmentation function D c ^y,{ z ) for c — > (cc)c — > (E cc , S* c )c is given by 
convolving the function D c ^( cc )(z) with the amplitude for the heavy antitriplet diquark to 
fragment to a ground state baryon with a single light quark. This latter function, which 
we shall denote D-q^ m (z), may be determined by data on charm and bottom antiquark 
fragmentation to heavy mesons. Summing over the E cc and S* c baryons, we then find 

A~e(*) = / —D c ^ (cc) (z/y)D^ M (y) . (3) 

j z y 

However, in the limit m c ^> Aqcd the heavy diquark carries all of the momentum of the 
baryon, and D-Q^ M (y) = 6(1 — y)Pq^ M , where Pq^ m is the integrated probability for a 
heavy antiquark to fragment to a ground state meson. Then (|3|) takes the simpler form 

D c ^(z)=Pq^ m D c ^ (cc) (z). (4) 
3 



Since in the limit we are considering all excited states will decay strongly to the ground 
state, we may replace this probability by unity. (There is a small correction from the 
fragmentation probability to baryons, Pq^x e> nere ^ s wou ld lead to an exotic ccqq final 
state.) Then the integral of (|l|) yields the final result 

/^ B „ EW = ± M 3™ c ) 2 ^^(H^-571n 2 ) . (5) 
Jo 9tt m% \ 30 J 

The E* c and E cc baryons will be produced in the ratio 2:1. However we may use the 
polarised fragmentation functions D c ^ cc \{z) derived in || to compute individually the 
populations of the various helicity states of the S cc and S* c . Let P be the net polarisation 
of the initial charm quarks. Then the populations of the baryon helicity states can be 
determined from P and the fraction £ of the produced diquarks which are transversely 
rather than longitudinally aligned. For example, a charm quark with energy E and helicity 
| can fragment to a (cc) diquark with helicity or 1, but not to one with helicity — 1, 
to leading order in m c /E. Hence the net polarisation of the diquarks is degraded to £P, 
where £ = 0.69 ||. To make the baryon, the diquark must then be combined with a 
light quark from the nonperturbative part of the fragmentation. In the limit m c ^> Aqcd, 
the helicity of the diquark is irrelevant to this soft process; the parity invariance of QCD 
then requires than the light quarks populate equally the helicities ±|. Since the E cc 
and Tj* c are produced incoherently, we may compute independently the probabilities that 
the light quark and the diquark will combine into the two possible angular momentum 
states, \ and |. We then find the various helicity states to be populated in the ratios 
E* c (±§) : E^ c (±±) : E cc (±|) = JC(1 ± P) : § + |C("1 ± |P) : |(1 ± (P). 

The quark fragmentation function (|3|) has been computed at the renormalisation scale 
fi = 3m c . The Altarelli-Parisi equations must then be used to evolve it up to a high scale 
\i = M typical of collider energies, a procedure which sums large logarithms of the form 
ln(m c /M). This evolution has two effects. First, it softens the z distributions, but because 
the relevant splitting functions P c ^ cg (z) and P c ^ gc (z) integrate to zero, the fragmentation 
probability L dz D c ^y,(z) remains unchanged. Second, a gluon fragmentation function 
D g ^-z(z, M) is induced via the gluon splitting function P g ^ C c(z). In leading logarithmic 
approximation, this effect, of order a^ln(M/m c ), dominates over any direct contribution 
to gluon fragmentation, which would be of order ct^ without the ln(M/m c ) enhancement. 

We now turn to the production of baryons of the form bcq, in which the heavy quarks 
are not identical. In this case the be diquark may be in either a spin-0 state, which we shall 
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denote (6c), or a spin-1 state, which we shall denote (be)* . In addition to the formation of 
Efo c and ££ c baryons from the fragmentation of the of the (be)* diquark, we now have the 
possibility of the formation of A& c baryons from the fragmentation of the (be) diquark. 

Just as in the case of diquarks consisting of two identical heavy quarks, it is straight- 
forward to relate the fragmentation functions for b — > (bc)^c and c — > (bc)^b to those for 
production of the physical states B c and B*. These latter functions have been calculated 
by Braaten, Cheung and Yuan |]], and we use their computation to obtain our results. Un- 
like the (cc) case, the colour wavefunction of the diquark is \(bc) a ) = -^e a bc \b b ) |c c ), and 
there is no factor of 2 from fermi statistics, so the spin-0 diquark fragmentation function 
with an initial bottom quark is given by 

D b ^ (bc) (z, fi ) = a s (no) 2 F(z, r) , (6) 

y 7T TTXj ^ 



where 



F(z, r) = 12(1 ^ (1 _ J r)z)6 [6 - 18(1 - 2r)z + (21 - 74r + 68r 2 )z 2 

-2(1 - r)(6 - 19r + 18r 2 )z 3 + 3(1 - r) 2 (l - 2r + 2r 2 )z 4 ] , 



r = m c /(m c + mj), and /j,q = ^Am c (mh + m c ) is a mass scale intermediate between m c 
and mj. (The selection of this scale is discussed in detail in 0.) The analogous function 
for an initial charm quark, -D c — >(6c)(' z ), is given by (|6|) with the replacements F(z,r) — > 
(m c /mf,) 3 F(z J 1— r) and /jlq — > /j' = sjAm^ra^ + m c ). Finally, the fragmentation function 
-Dfc^(bc)* (z) to the spin-1 diquark is given by © with F(z, r) — > F*(z, r), where 

F*(z, r) = " Z{ )~ Z) \ , 6 [2 - 2(3 - 2r)z + 3(3 - 2r + 4r 2 )z 2 

4(1 - (1 - r)z) b (8) 

-2(l-r)(4 - r + 2r 2 )^ 3 + (1 - r) 2 (3 - 2r + 2r 2 )^ 4 ] . 

Then the same replacements F*(z,r) — ► (m c /mb) 3 F*(z,l — r) and /io — ► A*o yield 
L> c ^ (6c) *(z). 

A calculation similar to that in Ref. gives £ = 0.69 for the net alignment of the 
(6c)*. Hence the populations of the various helicity states of the £fo c and ££ c will be in 
the same ratios as for the E cc and £* c . 

The diquark distributions ©-® must be convolved as in @ with experimentally 
determined meson fragmentation functions Dq^ m (z) to obtain fragmentation functions 
to Efe c , Y*l c and Aj, c - These distributions are then subject to Altarelli-Parisi evolution up 
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to collider energies. As in the case of (cc) production, the quark fragmentation functions 
are softened by this evolution, and gluon fragmentation functions are induced. The de- 
tailed effect of the Altarelli-Parisi equations on the fragmentation functions to B c and B* 
is presented in || , and the discussion given there applies here as well. Of course, the frag- 
mentation functions reduce directly to integrated probabilities as in (|J) and (^), quantities 
which do not evolve and are more accessible experimentally. 

It is interesting to note that at high-energy colliders the rates for production of these 
doubly heavy baryons are comparable to those for the more familiar quarkonium systems. 
Relating our fragmentation probabilities to those for c — > probability for c — > E cc , S* c 
to be ~ 2 x 10" 5 , for b -> A bc to be ~ 2 x 10" 5 , and for b -> £ &c , £* c to be ~ 3 x 10" 5 . The 
probabilities for b — > E^, ££ fe , c — > A& c and c — > E^c, E£ c are down by roughly (m c /mfo) 3 , 
or two orders of magnitude. Hence it may be possible to observe the doubly heavy baryons 

(*) A (*) 

£ cc j Afe c and E^/ at the Tevatron. 

There will be additional contributions to the fragmentation to doubly heavy baryons 
from radially excited diquark states which subsequently decay to the ground state. Equa- 
tions (U) and (H) also hold for these processes, using the appropriate value for the diquark 
wavefunction at the origin. Since for the ifi system, (R^(2S)/R^) 2 — 0.4, we expect that 
(R( cc (2S))/R(cc)) 2 is also n °t particularly small and that excited diquarks will contribute 
significantly to the production of doubly heavy baryons. 

The largest uncertainty in our calculation arises from our lack of knowledge of the 
diquark wavefunction at the origin. We have naively scaled the values for quarkonium 
systems by a colour factor of (1/2) 3 , which is valid only for wavefunctions living entirely 
in the Coulombic region of the potential. This is certainly a poor approximation for cc, 
bc and bb bound states. However, we note that in deriving (H) and (f|) we have assumed 
nothing about the potential except that the diquark is sufficiently tightly bound to have a 
~ 100% probability of fragmenting to baryons: f Q dy D {QQ) ^ QQq (y) = J Q dy Dq_^ m (y) = 
1. In reality, this integrated probability is somewhat less than one, since the diquark may 
dissociate as it hadronises, for example (cc) — > DD + X. Thus it is likely that (|5|) and © 
somewhat overestimate the true fragmentation probabilities. 

In summary, we have calculated the fragmentation functions of a heavy quark Q to a 
doubly heavy baryon of the form QQq or QQ'q. The perturbative part of the calculation 
can be related simply to the fragmentation function to quarkonium, and is of a similar 
magnitude. The nonperturbative part may be related to the fragmentation of heavy quarks 
to heavy mesons, which may be measured experimentally. We have used our results to 
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estimate fragmentation probabilities of b and c quarks to doubly heavy baryons of the form 
ccq, bcq and bbq. We find the production rates for doubly heavy baryons are large enough 
that some of these states may eventually be observable. 
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Fig. 1. Feynman diagrams responsible for the fragmentation c — > (cc)c. 
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This figure "fig2-l.png" is available in "png" format from: 
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